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Abstract 

We investigate the interplay between the symmetry and the emergence of the quarkyonic phase, 
adding the flavor-dependent complex chemical potentials fif = + ^TOf with {9f) = {0,0, —6) to the 
Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) model. When 6* = 0, the PNJL model with the ^ f 
agrees with the standard PNJL model with the real chemical potential /u. When 9 = 2tt/3, meanwhile, the 
PNJL model with the fi f has the Z symmetry exactly for any real /i, so that the quarkyonic phase exists 
at small T and large jj,. Once 9 varies from 27r/3, the quarkyonic phase exists only on a line of T = 
and fi larger than the dynamical quark mass, and the region at small T and large /i is dominated by the 
quarkyonic-like phase in which the Polyakov loop is small but finite. 
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Understanding of the confinement mechanism is one of the most important subjects in hadron 
physics. Lattice QCD (LQCD) shows numerically that QCD is in the confinement and chiral 
symmetry breaking phase at low temperature (T) and in the deconfinement and chiral symmetry 
restoration phase at high T. In the limit of infinite current quark mass, the Polyakov-loop is an 
exact order parameter for the deconfinement transition, since the Z^v^ symmetry is exact there. The 
chiral condensate is, meanwhile, an exact order parameter for the chiral restoration in the limit of 
zero current quark mass. In the real world where u and d quarks have small current masses, the 
chiral condensate is considered to be a good order parameter for the chiral restoration, but there is 
no guarantee that the Polyakov-loop is a good order parameter for the deconfinement transition. 

In the previous paper [jl|], we have proposed a QCD-like theory with the Z^r^ symmetry. Let us 
start with the SU(Ai'c) gauge theory with Nf degenerate flavors to construct the QCD-like theory. 
The partition function Z of the S\]{Nc) gauge theory is obtained in Euclidean space-time by 

Z = j DqDqDAexp[-So] (1) 

with the action 

^0 = 1 d'x[J2 QfiluD. + mj)qf + ^F^J, (2) 

where g/ is the quark field with flavor / and current quark mass nif, Dy = dy—iAy is the covariant 
derivative with the gauge field A^, g is the gauge coupling and F^^, = d^A^ — d^A^ — i[A^, Ay]. 
The temporal boundary condition for quark is 

g/(x,/3 = l/r) = -g/(x,0). (3) 

The fermion boundary condition is changed by the Zat^ transformation as [j2l |3|] 

/5) = - exp {-i2'Kk/ N^)qf{x, 0) (4) 

for integer k, while the action 5*0 keeps the form of (|2]) in virtue of the fact that the Zat^ symmetry 
is the center symmetry of the gauge symmetry [)2|]. The Z^r^ symmetry thus breaks down through 
the fermion boundary condition (|3]) in QCD. 

Now we consider the SU(A^) gauge theory with degenerate flavorSji.e. N = Nc = Nf, and 
assume the twisted boundary condition (TBC) in the temporal direction yj]: 



qf{x,P) = -exp{-i9f)qf{x,0) 
2 



(5) 



with the twisted angles 

% = 27r(/-l)/iV + ^i (6) 

for flavors / labeled by integers from 1 to A^, where 9i is an arbitrary real number in a range of 
< < 2n. The action Sq with the TBC is a QCD-like theory proposed in Ref. [1]. In fact, the 
QCD-like theory has the Z^r^ symmetry, since / is changed into / + A; by the Zat transformation but 
f + k can be relabeled by /. In the QCD-like theory, the Polyakov loop becomes an exact order 
parameter of the deconfinement transition. The QCD-like theory then becomes a quite useful 
theory to understand the confinement mechanism. 
When the fermion field qj is transformed by 

Qf — )■ exp {—i9fTT)qf (7) 

with the twisted angle 9f and the Euclidean time r, the action 5*0 is changed into 

8(9 f) = j d'x[J2 qAluD, - i9jT^, + mf)qf + ^F^J (8) 

with the imaginary chemical potential jJif = 'i'T9f, while the TBC returns to the standard one (|3j. 
The action So with the TBC is thus identical with the action S{9f) with the standard one (|3]). In 
the limit of T = 0, the action S{9f) comes back to the QCD action 5*0 with the standard boundary 
condition ([3]) kept. The QCD-like theory thus agrees with QCD at T = where the Polyakov 
loop is zero. One can then expect that in the QCD-like theory <P is zero up to some temperature 
Tc and becomes finite above Tc, i.e, that the Zat^ symmetry is exactly preserved below Tc but 
spontaneously broken above T^.. Actually, this behavior is confirmed by imposing the TBC on the 



211]. The PNJL model with the 



Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) modelyl- 
TBC [1] is referred to as the TBC model in this paper. In the TBC model, the flavor symmetry is 
explicitly broken by the flavor-dependent TBC dS]), but the flavor-symmetry breaking is recovered 
at T < Tc. The TBC model is thus a model proper to understand the confinement mechanism. 

A current topic related to the confinement is the quarkyonic phase lioi [ill [isl . 2o[ 221. It is 
a confined (color- singlet) phase with finite quark-number density n, that is, a phase with ^ = 
and n ^ 0. The n-generation induces the chiral restoration; in fact, the two phenomena occur 



almost simultaneously in the PNJL model 12111 . This fact indicates that the quarkyonic phase can 



be regarded as a chirally-symmetric and confined phase. It was suggested in Refs. [|23 



2411 that 



the chirally-broken phase is enlarged toward lager fi by the chiral density wave. In this paper. 



for simplicity, we ignore inhomogeneous condensates such as the chiral density wave. Effects 
of the inhomogeneous condensate on the quarkynic phase and the interplay between the effects 
and the TL^^ symmetry are interesting as a future work. The concept of the quarkyonic phase 
was constructed in large A^c QCD. In fact, the phase was first found at small T and large real 
quark-number chemical potential in large N^. QCD. Recently, the PNJL model showed that a 
quarkyonic-like phase with <P < 0.5 and n 7^ exists at small T and large pi, for the case of 



A^c = 3[113 



2011 . This result may stem from the fact that the deconfinement transition is crossover 
for Ai'c = 3. This suggests that the quarkyonic phase can survive even at Ai'c = 3 in the QCD-like 
theory with the TLjq^ symmetry. 



Fig. 1: Location of exp[z^?j] in the complex plane; here, (Qj) = (0, 6, —6). 

In this paper, we consider the PNJL model of = A^c = ^/ = 3 with the flavor-independent 
real chemical potential ji and the flavor-dependent quark boundary condition (|5} with 

{ef) = {Q,e,-e) (9) 

instead of Q; see Fig. [T]for the boundary condition. The present system is the same as that 
with the standard boundary condition (|3]) and the flavor-dependent complex chemical potentials 
jjLf = jjL + iTOf with (|9l). The present model with the yU/ is reduced to the standard PNJL model 
with the flavor-independent real chemical potential jj, when 6* = and to the TBC model with 
the ZiVc symmetry when 6 = 2n/3. Varying 9, one can see how the phase diagram is changed 
between the exact color-confinement in the TBC model and the approximate one in the standard 
PNJL model. The aim of this paper is to see this behavior. Our particular interest is the location 
of the quarkyonic and the quarkyonic-like phase in the fi-T plane. 

In general, there is no guarantee that the QCD partition function with complex chemical po- 
tential is real. It is, however, possible to prove that the QCD partition function Zo{nf) with 



fif = fi + iT9f is real. The fermion determinant det A^(/i/) satisfies the relation 



[detMif^f)]* = detM{-fx}) 

= Yl det[D -{fi- iOfT)-^4 + rrif] 
f 

= Yl det[D - (/i + z%T)74 + ruf] = det A^(-/i/) 
/ 



(10) 



where the third equality is obtained by the relabeling of the /. The present system thus has the 
sign problem, but the partition function is real, since 



(11) 



where the first equality is obtained by (flOl ) and the second one by the charge conjugation. Also in 
the PNJL model with the /i/, the partition funciton is real, as shown later. 
The three-flavor PNJL Lagrangian is defined in Euclidian space-time as 
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a=0 



detgi(l + 75)'?j + h.c. 



+ Ui<P[A],<P*[A],T), 



(12) 



where = d,^ — id^^A^, is the Gell-Mann matrices and m = diag(mi, m2, m-^). Gs and Gd 
are coupling constants of the scalar-type four-quark and the Kobayashi-Maskawa-'t Hooft (KMT) 



interaction n25[ 



2611 . respectively. The KMT interaction breaks the t/A(l) symmetry explicitly. The 



Polyakov-loop <P and its conjugate <P* are defined by 

^ = itre(L), ^* = ltre(L), 
with L = exp{iA4/T) in the Polyakov gauge. We take the Polyakov potential of Ref. ||8D 

a(r) 



(13) 



*\2\ 



-<P*<P + b{T) ln(l - 6$$* + 4(^^ + - 3{$$*) 



a{T) = ao + ai + aa ( ' ^(^) = ^3 (^) • 



(14) 
(15) 



Parameters of U are fitted to LQCD data at finite T in the pure gauge limit. The parameters except 
To are summarized in Table H The Polyakov potential yields the first-order deconfinement phase 
transition at T = Tq in the pure gauge theory [1271 l28ll . The original value of To is 270 MeV 
determined from the pure gauge LQCD data, but the PNJL model with this value yields a larger 
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ao 




02 


b3 


3.51 


-2.47 


15.2 


-1.75 



TABLE I: Summary of the parameter set in the Polyakov-potential sector determined in Ref. ||8|l. All 
parameters are dimensionless. 



value of the pseudocritical temperature Tc at zero chemical potential than Tc ~ 160 MeV predicted 
by full LQCD ibsUsil]. We then rescale Tq to 195 MeV so as to reproduce = 160 MeV jm. 

Now we consider the flavor-dependent complex chemical potential jjf = n + iO/T. The ther- 
modynamic potential (per volume) is obtained by the mean-field approximation as [|l6ll 



(16) 



with 



Q 



rf3 



P 



(27r)^ 



Ef + -^\n(l + Le-^'^f) + ^ In (l + Le"^^" 



(17) 



where o"/ = {qfQf), = Ef ± fif and Ef = y p"^ + Mf \ Here the three-dimensional cutoff is 
taken for the momentum integration in the vacuum term [|l6|l. Obviously, is real. The dynamical 
quark masses Mj and the mesonic potential Um are defined by 



Mf = mf - AGsCTf + 2Gd I e/gfe I cTgcr^i, 

= ^ ^Gso"/ - 4G'DO"io-2a3, 
/ 



(18) 
(19) 



where ejgh is the antisymmetric symbol. 

The PNJL model has six parameters, (Gs, Gd, mi, m2, m^. A). A typical set of the parameters 
is obtained in Ref. ll32ll for the 2-1-1 flavor system with mi = m2 = mi < m^. The parameter 
set is fitted to empirical values of r/'-meson mass and vr-meson mass and vr-meson decay constant 
at vacuum. In the present paper, we set m^ = mi = m^ in the parameter set of Ref. [|32|l. The 
parameter set thus determined is shown in Table im 

Taking the color summation in (fT6] ) leads to 



-2E 



P 



3Ef + - (In J^f + In J^j) J +UM{crf) + U{<P, T), 



(20) 
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mo(MeV) 


yl(MeV) 






5.5 


602.3 


1.835 


12.36 



TABLE II: Summary of the parameter set in the NJL sector. All the parameters except mo are the same as 



in Ref. 1321. 



where 

Tf=l + 3<Pe-'^^y + 3<P*e-^'^^f + e'^^^f , (21) 
J^f= l + 3<P*e-^^t +3(Pe-^^'^f +e-^^'^f . (22) 

Note that iJ^2) is the complex conjugate to (J-g), indicating that is real. 
In the case of 6* = 27r/3, particularly, H is invariant under the Z3 transformation, 

Namely, H possesses the Z3 symmetry. When the exact color-confinement with <P = occurs, f2 
is invariant for any interchange among E^, and Ef. Namely, (2 has the flavor symmetry in 
the exact color-confinement phase. 

Figure |2] shows T dependence of (a) the Polyakov loop ^ and (b) the chiral condensate ai at 
/i = 0. The solid, dashed and dotted curves represent three cases of = 0, Svr/lS and 27r/3, 
respectively. For 9 = corresponding to the standard boundary condition, the chiral and decon- 
finement transitions are both crossover. For 9 = 2n/3 corresponding to the TBC, the first-order 
deconfinement transition occurs at T = Tc = 203MeV and the exact color-confinement phase 
appears below Tc. The first-order transition of at T = Tc propagates to a/ as a discontinuity. 
For 9 ^ 2n/3, the deconfinement transition is no longer exact. As 9 decreases from 27r/3 to zero, 
T dependence of <P becomes slower, and near 9 = n/2 the order of the deconfinement transition 
is changed from the first-order to crossover. 

Figure |3] shows the phase diagram in the T-yU plane. Panels (a)-(c) correspond to three cases 
of 6* = 0, 87r/15 and 27r/3, respectively. The thick (thin) solid curve represents the first-order 
deconfinement (chiral) phase transition line, while the thick (thin) dashed curve corresponds to 
the deconfinement (chiral) crossover line defined by the peak of d^/dT (daj/dT). For 6^ = 0, 
the chiral and deconfinement transitions are both crossover at smaller /i, but the former becomes 
the first-order at larger /x. For 9 = 27r/3, the deconfinement transition is the first-order at any 
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Fig. 2: (a) The Polyakov loop ^ and (b) the chiral condensate a\ in the Q-T plane at /i = OMeV. 



/i, whereas the first-order chiral transition line appears only at ~ Mj = 323 MeV. The region 
labeled by "Qy" at /i > Af/ and small T is the quarkyonic phase, since <P = and n ^ there. 
The region labeled by "Had" is the hadron phase, because the chiral symmetry is broken there and 
thereby the equation of state is dominated by the pion gas plQ. The region labeled by "QGP" 
corresponds to the quark gluon plasma (QGP) phase, although the flavor symmetry is broken there 
by the TBC. As 9 decreases from 2n/3 to zero, the first-order chiral transition line declines toward 
smaller fi and the critical endpoint moves to smaller fi. Once 9 varies from 27r/3, the quarkyonic 
phase defined by ^ = and n ^ shrinks on a line with T = and /i > Mj and a region at small 
T and /i > A// becomes a quarkyonic-like phase with small but finite <P and n ^ 0; the latter 
region is labeled by "Qy-like". 

For small 9 far from 27r/3, the deconfinement transition line declines as ji increases, but for 
9 = 2n/3 the line is almost horizontal at small fi and rises at intermediate fi, as seen in Fig. |3l The 
rising of the deconfinement transition line is a consequence of the Z3 symmetry, as shown below. 
The quark one-loop part of H, which is defined by J7q in (flTI) . can be expanded into a Maclaurin 
series 



i7Q = i7Q((p = 0, <P* = 0)+ Cio^ + Coi<P* + C20^'' + Cu<P<P* + Co2^*^ + 



(24) 
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Fig. 3: Phase diagram in the T-/i plane. Panels (a)-(c) correspond to three cases of = 0, Stt/IS and 27r/3, 
respectively. The thick (thin) solid curve means the first-order deconfinement (chiral) phase transition line, 
while the thick (thin) dashed curve does the deconfinement (chiral) crossover line. The closed circles stand 
for the endpoints of the first-order deconfinement and chiral phase transition lines. In panels (a) and (b), the 
thick-solid line at T = and [i > Mj = 323 MeV represents the quarkyonic phase. 



The coefficients Cnm are explicitly obtained as 



ClO = - 



Coi 



< 0, 



< 0, 



C20 



Cll = -77 



C02 




> 0, 



> 0, 



> 0. 



(25) 
(26) 
(27) 
(28) 
(29) 



The Cnra STC positivc for even n + m but negative for odd n + m. The absolute values of the 
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Cnm increase as /i increases, unless fi is quite large. For simplicity, we fix Mf to a constant to 
focus our attention on $ dependence of H. In this assumption, Um and the zeroth-order term 
i7Q(<P = 0,<1>* = 0) in the Maclaurin series are just constants and thereby become irrelevant to 
the present discussion. So we neglect these terms. We also assume that (p = (p* . This is true for 
/i = and well satisfied for small and intermediate ji of our interest. In the pure gauge limit where 
i^Q = 0, the thermodynamic potential Q agrees with the Polyakov potential U{$) and hence has 
no ji dependence. The potential has a global minimum at <P = and a local one at = > 
for small T: namely, U{<P = 0) < U{<P = <Prn)- For the case of = 27r/3, the system has the Z3 
symmetry. Up to the second order of the Maclaurin series, only the cii<P<P* term appears because 
of the symmetry. When the term is added to hi, the resultant potential keeps the same value as 
U{$) at ^ = 0, but increases from U{(p) at <P > 0. This property makes the deconfinement 
transition more difficult. The coefficient cn as a function of fi little increases for fi <^ Mf, but the 
increase becomes sizable for /i > Mj. Therefore the rising of the deconfinement transition line 
with respect to increasing /i is tiny at small fi but becomes sizable at intermediate fi. The potential 
cii<P<P* + U keeps a positive curvature at <P = because of cn > 0, so that the deconfinement 
transition is first-order for any positive value of cn. For the case of ^ 2n/3, meanwhile, the 
first-order term cio<? + Coi<P* is not prohibited by the Z3 symmetry and thereby dominates Hq 
particularly for small 9 far from 2n/3. Since cio and cqi are negative, the situation for small 9 
becomes opposite to that for 9 = 2n/3. Eventually, the deconfinement transition line slopes down 
as /i increases for the case of small 9. 

Figure |4] shows T dependence of <1> for /i = 0.1, 0.3 GeV. We consider the case of 6* = in 
panel (a) and that of = 2n/3 in panel (b) by assuming f] = cio<? + cqi^* + U in panel (a) 
and H = cii<p(p* + U m panel (b). Note that Mf is fixed to 323 MeV and T dependence of <P is 
determined from the Q with the minimum condition. As fi increases, the transition temperature 
decreases for = 0, but increases for 9 = 2n/3. The transition is first-order for the case of 
9 = 271 / 3. These results are consistent with the qualitative discussion mentioned above. 

In summary, we have investigated the interplay between the Ztv^ symmetry and the emergence 
of the quarkyonic phase, adding the complex chemical potentials fif = jj, + iT9f with {9f) = 
(0, 9, -9) to the PNJL model. When 9 = Q, the PNJL model with the yU/ is reduced to the PNJL 
model with real /x. This situation corresponds to QCD at real fi. When 9 = 2n/3, meanwhile, 
the PNJL model with the yU/ is reduced to the TBC model with the Z^^ symmetry. This situation 
corresponds to the QCD-like theory with the Z^v^ symmetry at real jj,. When 9 = 2n/3, the 
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T (GeV) T (GeV) 

Fig. 4: T dependence of ^ for = 0.1, 0.3 GeV in tlie lowest order approximation. Panel (a) corresponds 
to the case of ^ = and panel (b) does to the case ofO = 2n/3. 

quarkyonic phase defined hy <P = and n > really exists at small T and large /i. Once 9 
varies from 27r/3 to zero, the Zjv^ symmetry is broken. As a consequence of this property, the 
quarkyonic phase exists only on a line of T = and fx > Mf, and the region at small T and large 
fi is dominated by the quarkyonic-like phase characterized by small but finite <P and n > 0. The 
Ztv^ symmetry thus plays an essential role on the emergence and the location of the quarkyonic 
phase in the fi-T plane, and the quarkyonic-like phase at 6* = is a remnant of the quarkyonic 
phase at = 27r/3. Since the Z^^ symmetry is explicitly broken at = 0, it is then natural to 
expand the concept of the quarkyonic phase and redefine it by a phase with small ^ and finite 
n. For this reason, the quarkyonic-like phase is often called the quarkyonic phase. The gross 
structure of the phase diagram thus has no qualitative difference between 6 = 2tt/3 and zero, if 
the concept of the quarkyonic phase is properly expanded. In this sense, the Zjy^ symmetry is a 
good approximate concept for the case of 6* = 0, even if the current quark mass is small. 
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